Deriving local order parameters from tensor network representations 
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A large class of quantum phase transitions for quantum lattice systems are characterized by local 
order parameters. It is shown that local order parameters may be systematically constructed from 
tensor network representations of quantum many-body ground state wave functions by investigating 
the reduced density matrices for local areas on an infinite-size lattice. Depending on whether or not 
the system symmetries are spontaneously broken, and whether or not the ground state fidelity per 
lattice site is continuous, there are four categories of quantum phase transitions for systems with 
local order parameters. Quantum phase transitions characterized by nonlocal order parameters are 
discussed, aiming at better understanding quantum systems with topological order. 
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Quantum phase transitions (QPTs) [H are driven by 
quantum fluctuations due to the uncertainty principle, 
and occur at zero temperature when some control pa- 
rameter varies. In the last decades, much attention has 
been paid to QPTs in strongly correlated electronic sys- 
tems, since some of the most exciting discoveries in con- 
densed matter physics, such as various magnetic order- 
ings, the integer and fractional quantum Hall effects, 
and the high-T c superconductors, are often attributed 
to quantum critical phenomena. In the conventional 
Landau-Ginzburg- Wilson paradigm, the key notions are 
symmetry-breaking orders and local order parameters, 
with symmetry-broken phases characterized by the non- 
zero values of local order parameters. A new paradigm 
emerges for systems with topological order where no 
local order parameter exists to characterize exotic orders. 

In spite of their decisive role in characterizing QPTs, 
no scheme, which is applicable (at least in principle) 
to quantum systems undergoing QPTs, is available to 
systematically construct order parameters (either local 
or nonlocal). The only exception, to the best of our 
knowledge, is the work by Furukawa, Misguich and Os- 
hikawa |3| , who attempted to systematically derive order 
parameters from comparing the reduced density matri- 
ces of the degenerate ground states for various subareas 
of the system defined on a finite-size lattice. In their ap- 
proach, the low-lying energy spectrum and eigenvectors 
of a finite-size system with a given discrete symmetry are 
known (analytically or numerically), where a finite num- 
ber of ground states are assumed to be nearly degenerate, 
with their quantum numbers indicating what symmetries 
are broken in the thermodynamic limit. Therefore, the 
approach, as it stands, only works for systems with a 
finite number of (nearly) degenerate ground states. 

The difficulties of determining order parameters for a 
QPT in a given quantum system lie in the facts that 
(i) usually it is a formidable task to compute ground 
state wave functions and (ii) the ground state phase dia- 
gram is often lacking. However, significant advances have 



been made in both the classical simulation of quantum 
lattice systems and the determination of ground state 
phase diagrams: first, a tensor network (TN) represen- 
tation of quantum many-body wave functions provides 
an efficient way to classically simulate quantum many- 
body systems 0, H @, 0, H, IE EE El ; second, two novel 
approaches to study QPTs have been proposed from a 
quantum information perspective, namely entanglement 
13, EE 0M [IE El El and fidelity [IE SE III IS E3 • 
In Ref. [22|, a viable scheme to determine the ground 
state phase diagram of a quantum lattice system with- 
out prior knowledge of order parameters was proposed. 
This was achieved by studying the singularities in the 
ground state fidelity per lattice site [13| , combining with 
a practical way to compute the fidelity per lattice site for 
infinite-size lattice systems using the TN algorithms for 
translationally invariant systems [E Q . 

In this paper, we develop a practical scheme to system- 
atically determine local order parameters for quantum 
many-body systems on infinite-size lattices from tensor 
network representations of the ground state wave func- 
tions by investigating the reduced density matrices for 
local areas. It yields that there are four categories of 
QPTs characterized by local order parameters, depend- 
ing on whether or not the ground state fidelity per lat- 
tice site is continuous, and whether or not the system 
symmetries are spontaneously broken. If no local order 
parameters exist, then the system is described by a non- 
local order parameter. The latter is relevant to QPTs in 
quantum systems with topological order. 

Local and nonlocal order parameters. Consider a trans- 
lationally invariant infinite-size quantum lattice system 
S in D spatial dimensions described by a Hamiltonian 
H(X) with a global symmetry group G, where A is a 
control parameter Suppose the system undergos a 

QPT at a transition point A = A c . Then any two differ- 
ent ground states |^(Ai)} and l^^)} corresponding to 
two different values Ai and A2 of the control parameter 
A are orthogonal to each other, i.e., (ip(X2)\ip(Xi)) = 0. 
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Therefore, the ground states for different values of 

the control parameter A are reliably distinguishable via 
quantum measurements 25] . This leads to the existence 
of local physical observables, which enable to distinguish 
the states by means of local measurements, if sufficient 
copies of the system are simultaneously prepared [26| for 
different values of the control parameter. However, a 
more pertinent question is to ask if there is any local 
physical observable O (Hermitian and traceless) defined 
in a local area Q on the lattice that may tell us in which 
phase a given state |^(A)) is: (O) ^ in one phase and 

(0) = in the other phase, if copies of the system for one 
single value of the control parameter are prepared [2^ |. 
If the answer to the question is affirmative, then we may 
define the local physical observable as a local order pa- 
rameter [28} . Depending on whether or not the fidelity 
per lattice site is continuous [2^ |. and whether or not 
the symmetries are spontaneously broken, we have four 
different categories of QPTs with local order parameters: 

(1) discontinuous, symmetries not broken; (2) discontinu- 
ous, symmetries broken; (3) continuous, symmetries not 
broken; and (4) continuous, symmetries broken. If no 
local order parameter exists, then one may further seek 
a nonlocal physical observable O n defined on a nonlocal 
area on the infinite-size lattice to judge in which phase 
a given state \ip(X)) is, with only copies of the system 
prepared for one single value of the control parameter. 
Such a nonlocal physical observable O n is a nonlocal or- 
der parameter. This accommodates systems with exotic 
topological order [3(| , which may even coexist with local 
symmetry-breaking orders. 

Reduced density matrices for local areas on an infinite 
lattice. Now we are in a position to clarify what conse- 
quences one may draw from the existence of a local order 
parameter (O) on a local area fl. Here we emphasize 
that, due to translational invariance, what really mat- 
ters is only the size and shape of the area Q, rather than 
its relative position in the entire lattice. For our pur- 
pose, we partition the whole system into two parts-the 
local area Q and its complement to the entire lattice f2. 
As such, the local area is described by the reduced den- 
sity matrix pn(X) corresponding to a given ground state 
(■0(A)). Generically, the reduced density matrix pn(A) ex- 
hibits nontrivial form due to the distinguishability by lo- 
cal measurements 26[ , except for the constraints imposed 
by the global symmetries (in the symmetric phase). An 
important fact is that the reduced density matrix pn(A) 
possesses different structures in two phases as far as the 
nonzero entries are concerned: po(A) = poo (A) in one 
phase and p n (X) = poo(A) + (0)O/(TrO 2 ) in the other- 
phase, as follows from the presence of the local order pa- 
rameter (O). Two cases should be distinguished: (i) if 
O is not invariant under the symmetry transformation, 
then the symmetry is spontaneously broken; (ii) if O is 
invariant under the symmetry transformation, then no 
symmetry is broken. 



We stress that local order parameters are not unique. 
Indeed, suppose a local order parameter (O) exists for a 
local area O, then there exists another local order param- 
eter (O) for a larger area Cl D f2, since all information 
encoded in the reduced density matrix pn(X) should be 
encoded in the reduced density matrix Pq(X), due to the 
fact that pn(X) is obtained from Pq(X) by tracing out ex- 
tra degrees of freedom in the complement set f2/0. This 
implies that there is an optimal local order parameter 
that corresponds to the smallest local area for the sys- 
tem considered [3l[. 

On the other hand, if the reduced density matrices for 
all possible local areas share the same nonzero entries 
structure in two phases, then a nonlocal order parameter 
is necessary to characterize the QPT. That is, an area 
with nontrivial topology should be chosen to see if the 
corresponding reduced density matrices exhibit different 
structures in different phases [32| . 

Computation of the reduced density matrices from ten- 
sor network representations. Now we show that it is fea- 
sible to compute the reduced density matrices for vari- 
ous local areas on infinite-size lattices. In this regard, 
we rely heavily on the fact that the TN representation 
of quantum many-body wave functions provides a pow- 
erful means to efficiently simulate infinite-size quantum 
lattice systems in one [8| and two and higher |9( spatial 
dimensions. In one spatial dimension, the NT represen- 
tation is the matrix product states (MPS) 0], and, in two 
and higher spatial dimensions, the NT representation is 
the projected entangled-pair states (PEPS) Q. In both 
cases, the TN algorithms for infinite-size lattices offer an 
efficient way to compute the reduced density matrices 
for various local areas, if quantum lattice systems are in 
gapful phases. 

If quantum lattice systems are in gapless phases, then 
a more sophisticated representation, i.e., the multi-scale 
entanglement renormalization ansatz (MERA) fiol . 11 | is 
needed for quantum many-body states. A MERA also 
provides an efficient representation of quantum ground 
states on an infinite-size lattice. Thanks to the fact that 
the width of the causal cone is bounded, a MERA offers 
an efficient way to compute the reduced density matrices 
for various local areas. In Ref. it is described explic- 
itly how to compute the one-site and two-site reduced 
density matrices from a MERA. 

A practical scheme to derive local order parameters 
from tensor network representations. For a quantum lat- 
tice system with a symmetry group G, one may system- 
atically derive local order parameters from the TN rep- 
resentations of quantum many body ground state wave 
functions. It consists of two steps: (1) determine the 
ground state phase diagram by computing the ground 
state fidelit y p er lattice site in terms of the infinite TN 
algorithms [22]; (2) derive local order parameters from 
a representative ground state wave functions by investi- 
gating the reduced density matrices for local areas on an 
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infinite-size lattice. 

To be self-contained, let us briefly recall the defini- 
tion of the ground state fidelity per lattice site d(Xi, A 2 ). 
For any two ground states and |V->(A 2 )), the fi- 

delity F(Ai,A 2 ) asymptotically scales as i^(Ai,A 2 ) ~ 
d(^i,\2) N , w hh N the total number of sites in the lat- 
tice. Here d(Ai, A 2 ) is the fidelity per lattice site [20l.l22l|. 
which is well defined in the thermodynamic limit: 



d(A l7 A 2 ) = lim F-(A!,A 2 ). 

N— too 



(i) 



It satisfies the properties inherited from the fidelity 
^(Ai,A 2 ): (i) normalization d(Ai,Ai) = 1; (ii) symme- 
try d(Xi, A 2 ) = d(A 2 ,Ai); and (iii) range < d(Ai,A 2 ) < 
1. In fact, the ground state fidelity F(Ai,A 2 ) may be 
mapped onto the partition function of a £>-dimensional 
classical statistical vertex lattice model with the same ge- 
ometry [22] . Thus, the fidelity per lattice site rf(Ai, A 2 ) is 
nothing but the partition function per site in the classi- 
cal statistical vertex lattice model [33j . This justifies why 
QPTs may be detected as singularities in lnd(Ai, A 2 ) as 
a function of Ai and A 2 . 

Once the ground state phase diagram is determined, 
one may choose a representative ground state from each 
phase and investigate the reduced density matrix for a 
local area fi on an infinite-size lattice. If the nonzero 
entries structures of the reduced density matrices for a 
given local area f2 are different for different phases, then 
one may read off a local order parameter: (i) if the re- 
duced density matrices are invariant under the symme- 
try group G, then no symmetry is spontaneously broken 
(and the ground state is non-degenerate); (ii) if one of the 
reduced density matrices is not invariant under the sym- 
metry group G, then the symmetry is spontaneously bro- 
ken (and degenerate ground states arise). If the nonzero 
entries structure of the reduced density matrices for all 
possible local areas f2 is the same for different phases, 
then there is no local order parameter (34| , and vice versa. 
For systems with topological order, this is consistent with 
the fact that different states in the ground state subspace 
share the same bulk tensors in the TN representations, 
with the only difference at the top tensor [35| . 

Examples. Let us give a few examples to illustrate our 
general scheme. The first example is the two-dimensional 
quantum Ising model described by the Hamiltonian: 



(2) 



(r,r') r r 

Here a x and a z are the Pauli matrices at the lattice 
site r, and the control parameters A and e correspond 
to transverse and parallel magnetic fields. Note that if 
e = 0, then the model enjoys the Z 2 symmetry; oth- 



erwise no symmetry exists. As shown in Ref. [22|, for 
e = 0, the fidelity per site d(Ai,A 2 ) is continuous, but 
it exhibits a pinch point singularity, indicating a con- 
tinuous phase transition at A c « 3.044 36]. On the 



other hand, the one-site reduced density matrix p\ dis- 
plays different nonzero-entries structures in two phases: 



['1 
Pi 



l/2{a l f l )a z , with (aP) being 



1/2 + l/2{ai')a 
zero for A > A c and nonzero for A < A c , as evaluated 
using the infinite TN algorithm in Ref. [j|. p\ is non- 
trivial under the global Z 2 symmetry transformation, so 
the Zi symmetry is spontaneously broken. This implies 
the existence of a local order parameter: (O) = (o~P). 
For nonzero e, when A < A c , the fidelity per site d(ei, e 2 ) 
is discontinuous [22j], implying that a first order QPT 
occurs when e changes sign. The one-site reduced den- 
sity matrix p\ displays different structures in different 
phases: (o, ) changes sign, since the symmetry group is 
trivial [23]. 

The second example is a spin 1/2 model with three- 
body interactions: 

H = £ 2(. 9 2 -l)« +1 -(l+. 9 ) 2 < + ( 5 ~l)V.f< +1 a.f +2 . 

i 

(3) 

It is ^-symmetric under the global spin reversal in the 
z direction. As emphasized in Ref. [2fj| . the param- 
eter space should be compactified by identifying g = 

+oo with g = —oo, due to the fact that H(+oo) = 

H(—oo). Since the ground state is an MPS [371 ] . 
one may extract the fidelity per site d as d(g,g') — 
V 1 + \99'\l\/i l + MX 1 + iFD ir 9 and 9' are in different 
phases, and d(g,g') = (1 + yW)/ Vi 1 + MX 1 + W\) 
if g and g' are in the same phase. Thus there are two 
transition points: g — and oo. All states for positive g 
flow to the product state (<? = 1) with all spins aligning 
in the x direction, and all states for negative g flow to 
the cluster state [38| (g — —1). Since the ground state is 
unique, so no symmetry is spontaneously broken. There- 
fore, the reduced density matrices should be invariant 
under the Zi symmetry group. Actually, the one-site re- 
duced density matrix p\ exhibits different nonzero entries 
structures in different phases: pi(i) = 1/2 + 1/2 (<rf )a x , 
with (erf) = Ag/(l + g) 2 for g > and (erf) = for g < 0, 
thus leading to the local order parameter (erf). 

The third example is the spin 1 XXZ model with uni- 
axial single-ion anisotropy described by the Hamiltonian: 



H = zZ 



,°i D i+1 



JzS; s, 



i+lj 



(4) 

where Sf (a — x, y, z) are the spin 1 operators at the 
lattice site i, and J z and D are the Ising-like and single- 
ion anisotropics, respectively. The model exhibits rich 
phases [39, 40, [4l|, |42|. For J z > 0, three gapful phases, 
i.e., the large-D, the Haldane, and the Neel phases oc- 
cur, which have the symmetric, fully broken and par- 
tially broken Z 2 x Z 2 symmetry, respectively. It has 
been argued that there is a tricritical point (Jt,D t ) ~ 
(3.20, 2.90) jUSi]. H Jz < Jt, then there are two critical 
values D Cl and D C2 when D varies from — oo to oo, char- 
acterizing the Ising-like transition from the Neel phase to 
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the Haldane phase and the Gaussian transition from the 
Haldane phase to the large-Z? phase, respectively. Be- 
yond the tricritical point where the Haldane phase dis- 
appears, the large- D-Neel transition was believed to be 
first order, although no final proof is available 42]. This 
has been confirmed numerically by evaluating the fidelity 
per site based on the infinite TN algorithm 43] . The one- 
site reduced density matrix p\ exhibits different nonzero- 
entries structures in different phases: it is invariant under 
the spin reversal Zi symmetry in the large- D phase, but 
not in the Neel phase. This indicates that the Zi symme- 
try is spontaneously broken, with a local order parameter 
(O) = (—lY(Sf). No local order parameter exists in the 
Haldane phase [441 ] . 

Support from the Natural Science Foundation of China 
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